Time-dependent barrier passage of Two-dimensional non-Ohmic damping system 
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The time-dependent barrier passage of an anomalous damping system is studied via the general- 
ized Langevin equation (GLE) with non-Ohmic memory damping friction tensor and corresponding 
thermal colored noise tensor describing a particle passing over the saddle point of a two-dimensional 
quadratic potential energy surface. The time-dependent passing probability and transmission coef- 
ficient are analytically obtained by using of the reactive flux method. The long memory aspect of 
friction is revealed to originate a non-monotonic (5(power exponent of the friction) dependence of 
the passing probability, the optimal incident angle of the particle and the steady anomalous trans- 
mission coefficient. In the long time limit a bigger steady transmission coefficient is obtained which 
means less barrier recrossing than the one-dimensional case. 
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I. INTRODUCTION 



The problem of escape from a mctastable states poten- 
tial is ubiquitous in almost all scientific areas. A great 
amount of chemical events, such as chemical reactions, 
molecular diffusion, or collision of molecular systems, etc, 
can be modeled by a single barrier escape process within 
the framework of standard Brownian motion [J, [3) 9- 
Although many other thoughts such as transition state 
theory [1, d, ll] and unimolecular rate theory 0, H, Q , 
etc, still work, the celebrated landmark elucidation of 
this problem is the Kramers rate theory p^. Where in 
his famous work, H. A. Kramers established a reaction 
rate formula which is applicable to all the cases from 
moderate to strong damping. The method he used to 
calculate the rate constant is conventional, namely, flux- 
over-population ll| . However, we noticed, am ong the 
various theoretical concepts for rate calculation [ill [l^ , 
the method of reactive flux 13, 14, 15] is a most power- 
ful and convenient way to follow. In its spirit, the initial 
conditions are assumed to be at the top of the barrier, 
which correspond to the ensemble of trajectories which 
start with identical initial conditions but experience dif- 
ferent stochastic histories. The escape rate then can be 
calculated by investigating the various flux of particles 
passing the transition state at the saddle point. Achieve- 
ments of its application in one-dimensional (ID) cases 

[HQ 

ITU has been witnessed in the past decades, but 
it has never been used in the study of higher dimensional 
systems. 



The recent widespread interest in anomalous diffusion 
[H, [13, [3 has highlighted the critical role that it plays 
in characterizing a large group of nonstandard statistical 
physics. For which, the non-Ohmic model with a idiosyn- 
crat ic p ower law frequency depending spectral density 
[igl [20I [2H bears prominent consequences in describing 
a rich variety of frequency-dependent damping mecha- 



nism. Previous ID studies on it have shown that in the 
non-Ohmic case, the friction of the system has a long 
time memory, the transmission coefficient of the system 
acts as a non-monotonic function of (5(power exponent of 
the friction) which reveals a relatively strong recrossing 
phenomenon [H, [l^ . But no higher dimensional study 
has occurred on this subject. However, a recent work by 
us has shown that in high dimensional systems, the non- 
diagonal correlation between various degrees of freedom 
will endow the barrier passage a riveting character [23 |. 
in particular, particles moving in a two-dimensional (2D) 
potential energy surface (PES) tend to select an optimal 
diffusion path to surmount the barrier. So it is of great 
interest to investigate the time-dependent barrier passage 
of high dimensional non-Ohmic damping systems. 

The primary purpose of this paper is then to report our 
recent study on the non-Ohmic damping system. In Sec. 
II, the analytical expression of the saddle-point passing 
probability is obtained by solving the 2D coupled gener- 
alized Langevin equation with non-Ohimic friction ten- 
sor. In Sec. Ill, we give the time-dependent transmission 
coefficient derived by using of the reaction flux method. 
Sec. IV serves as a summary of our conclusion. 



II. 2D NON-OHMIC PASSING PROBABILITY 



We consider the directional diffusion of a particle in the 



2^1] 



with 



1,2 



2D quadratic PES: U{xi,X2) 

and detojy < 0, the motion of the particle is described 
by the generalized Langevin equation (GLE) with a non- 
Ohmic memory friction tensor: 



mijXj{t)+ dt'f3i-j{t-t')xj{t')+ujtjXj{t) =^ ^i{t), (1) 
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where the Einstein summation convention is used and 
the components of the random force are assumed to be 
zero-mean and their correlations obey the fluctuation- 
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dissipation theorem 



2KT r^J^ 



cosuj{t~t'), (2) 



with kg the Boltzmann constant, T the temperature of 
the reservoir and Jij{uj) takes the form of non-Ohmic 
spectral density Jij(w) — /uJrY ^ where 5 is the 

power exponent taking values between and 2, 7^- is the 
friction constant, and ujr denotes a reference frequency 
allowing for the damping constant 7ij to have the di- 
mension of a viscosity at any 5. 

The initial condition of the particle is denoted as 
Xj{Q) = XjQ and ij(0) = VjQ, where xiq < and 
vio > 0. Assuming that xi-axis is the transport direc- 
tion [cjii < 0], the reduced distribution function of the 
particle for xi, when the variables X2{t), vi{t) and V2{t) 
are integrated out, can be written as: 



, t, a^^^ , x^Q , v^Q , v^Q ) 

1 / ix^-{x^it))r 



2-K(Jj:^ [t) 



exp 



(3) 



Applying on it the Laplace transform technique to Eq. ([T|) 
we get xi{t) and its variance cr^j(t) at any time 



2 „t 

xi{i) = / H,{t - t%,{t')dt\ (4) 

<W - / dtiH,{t-ti) [ dt2Mtl)ij{t2))Hj{t~t2), 

Jo Jo 

(5) 

where the mean position of the particle along the trans- 
port direction is given by 



(6) 



which relates to the initial position and velocity. The 
time-dependent factors in Eq.® according to the resid- 
ual theorem are C^{t) = £-^[F,{s)/P{s)] {i = 1...4) 
with exponential forms, the two response functions in 
Eqs.Q and jS]) read Hi{t) = C-^[F5{s)/P{s)] and 
H2{t) = £~^[Fq{s)/P{s)], where denotes the inverse 
Laplace transform. The expressions of P{s) and Fi{s) 
(i — 1, ...,6) are written in the appendix. 

The probability of passing over the saddle point [xi — 
X2 = 0] is always called the characteristic function 
xix-^gjX^gjVj^gjV^glt) w^ch cau be determined mathe- 
matically by integrating Eq.Q over xi from zero to in- 
finity 



W{xi,t]x^„,x^^,v^^,v^^)dxi 

(7) 
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FIG. 1: Stationary passing probability as a function of the 
incident angel for various 5. The parameters used are uau = 
—0.5, mi2 = 0.6, LUii = —2.0, 1022 = 1.5, mn = 1.5, 77122 = 
2.0, 712 — 1.0, 711 — 1.8, 722 = 1.2, uir — 1.0, and the initial 
condition of the particle is set x\o — —1.0, a;2o = 0, t;o = 1.9. 



In the barrier surmounting process, the probability den- 
sity function of the particle that is originally at one side 
of the barrier spreads in the phase space due to thermal 
fluctuation, and its center moves driven by its initial ve- 
locity. For large times, the probability of passing over the 
saddle point converges to a finite value depending on the 
initial position and velocity of the particle, and the other 
fraction of the particles goes to infinite along the oppo- 
site direction. This results in the characteristic function 
Eq. ([7]) to be a number which is equal to 1 for reactive tra- 
jectories and for nonreactive trajectories. If we consider 
only the long time result of state transition induced by 
the particle diffusion, the stationary passing probability 
is of great important to rely on. Mathematically it is cal- 
culated by Ppass = limt^oo |erfc [-{xi{t)) /{V2ax,{t))]. 

In Fig. [U we plot the stationary passing probability as 
a function of the incident angel for various power expo- 
nents S. In the calculations here and following, we rescale 
all the variables so that the dimensionless unit such as 
ks — 1.0 is used. Each friction strength jij is fixed to be 
a constant independent of S, and the reference frequency 
uJr = 1-0. It is seen from Fig. [T]that, the stationary pass- 
ing probability of non-Ohmic case {6 = 0.6 or S = 1.8) 
is not so large as the result of Ohmic case {S — 1.0). 
This can be understood from the point of view of the 
particle's critical velocity which corresponds to the con- 
dition: limt^oo(a^i(^)) = 0. Supposing v-^^g = VgCOS(j) and 
v^g — VgSm(j) are the two components of the particle's 
velocity along two directions xi and X2- The critical ve- 
locity for a particle reads: 



^ Ci{co)x,„ +C2{oo)x^„ 
° C3 (oo)cos0 -I- C4 (oo)sin(/) ' 

For the parameters in Figll] we obtain 
2.6417 > v'ls^o.e = 2.3477 > v'Js^i.o = 



(8) 



1.8426 when 
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the incident angle is set (f> = 0.258rad. Thus we can see 
the particle of non-Ohmic damping system has a small 
passing probability. 

It is also shown in Fig. [T]that, for each 5 there exits 
a certain incident angle for the particle to reach a max- 
imum stationary passing probability. This implies, for 
the 2D non-Ohmic damping system, an optimal incident 
angle (or an optimal diffusion path) still exists for the 
particle to surmount the PES, in accordance with the 
Ohmic case {5 ^ 1.0) we have studied [22|. The optimal 
incident angle for a particle of 2D non-Ohmic damping 
system can also be determined by retrospect the mini- 
mum value of the critical velocity Eq. ([5]). After some 
algebra we find in the long time limit it can be expressed 
by the system parameters as 
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)FIG. 2: Optimal incident angel as a function of the power 
iQ\ exponents 5 for various effective frictions. The different pa- 
rameters for each curve are (a) 712 = 0.8, 711 = 1.8, 722 = 1-2; 
(b) 712 = 2.0, 711 = 2.5, 722 = 2.2; (c) 712 = 4.0, 711 = 4.8, 
722 = 4.2, others are the same as those used in Fig. 1. 



where e is the largest analytical root of P(s) = 0, 
~ is the Laplace transform of the friction 

kernel function and 7ij = ^ijijj}.~^smr'^(5i: /2) is namely 
the effective friction here. 

Noticing that Eq. ([9]) implicitly implies the power ex- 
ponent 5. This means in the non-Ohmic damping case 
the optimal incident angle has an intimate relation to 
the friction of the system. In Fig. [2l the 5 dependence of 
the optimal incident angle ip at various effective friction 
is plotted. It shows that ip evolves as a non-monotonic 
function of 5. This implies that the barrier passage of 
the two-dimensional non-Ohmic damping system is inti- 
mately controlled by the friction. Particle facing different 
friction strength will select different optimal path in its 
barrier surmounting process. 

For instance Fig. [3] gives a schematic illustration of 
the optimal incident angel for various power exponents 
5 at different friction strength. From which we can see, 
comparing with the Ohmic case (i.e. 5 = 1.0), when the 
system friction strength is not so strong (straight lines 
in Fig. [3]), the optimal incident angle of non-Ohmic case 
tends to approach zero from a clockwise direction at most 
values of 5. However, if the system has a relatively strong 
friction (dashed lines in Fig. [3]), the optimal incident an- 
gle of non-Ohmic case tends to approach zero from a 
anticlockwise direction. Since the direction of zero inci- 
dent angle represents the valley direction of the PES, this 
implies the particle diffusion under the influence of non- 
Ohmic friction will has the probability of tracing along 
the valley passage. This is a meaningful prediction for 
the study of barrier escape problem. To deicide the re- 
lation between the system parameters that enables the 
particle to travel along the valley direction, one needs 
only set -0 = 0. This will be very useful in studying the 
fusion reactions of massive nuclei. 
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FIG. 3: Schematic illustration of the optimal incident angel 
for various power exponents 5. The different parameters for 
each group of curves are (a) 712 = 0.8, 711 — 1.8, 722 = 1.2 
for straight lines; and (b) 712 = 4.0, 711 = 4.8, 722 = 4.2 for 
dashed lines, other parameters are the same as those used in 
Fig. 1. 



III. NON-OHMIC ANOMALOUS 
TRANSMISSION COEFFICIENT 

In this section we evaluate the rate of a particle escape 
from the 2D metastable potential. Since difficulty makes 
it impossible to solve the Fokker-Planck equation with 
non-Ohmic friction, we use the method of reactive flux 
to derive the rate constant from Langevin dynamics. In 
the spirit of reactive flux calculation, particles are set 
evolving from the top of the barrier and the thermal rate 
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constant is got from the ensemble average of trajectories 
starting with identical initial conditions but experience 
different stochastic histories. For a 2D system, the rate 
for a particle of unit mass [1^ is defined as 

m = ^/ dpas(p.n,)x(p,g;)e-^(^^''-=)/'=B^. (10) 

Where p is the momentum vector and Q is the partition 
function for reactants integrating over the distribution 
of the ground states. x(Pj9s)) namely the characteristic 
function is a number which is equal to 1 for reactive tra- 
jectories {p-fis > 0) and for nonreactive ones {p-fis < 0). 
Because of the stochastic nature of the dynamics, it is 



necessary to take into account all the different possible 
realizations of a trajectory. To do so x(p, (js) should be re- 
placed by its non-equilibrium average xi^ia 1 2^20 1 '^lo 1 ^20 ) 
which results in a weighted factor with its value rang- 
ing from to 1. The rate is then proportional to the 
total flux from reactants to products independent of the 
choice of dividing surface between reactants and products 
excluding the influence of the recrossing effects. 

If the initial conditions are assumed to be at the top of 
the barrier, the rate constant corresponds to an ensemble 
of particles starting from (x^^, =^ 0, x^^, v-^^jV^^ ) at t = 
and obeys the equilibrium distribution. Mathematically, 



1 ^ f°° 



i=l 
'V„/kBT 2 

n 

i=l 



Qh 



(11) 



^rriijV.gV^g + ^ujijX.^x.Q + Vg , and an initial Boltzmann form 



with the system Hamiltonian iJ(a;j(j, x^^, , Ujo) 
stationary probability distribution with a weight ksT at the metastable well is assumed when the temperature is 
much less than the barrier height of the metastable potential. 

The total rate thus can be viewed as a transition state theory (TST) rate fc^-^^ = ^s"^^^'^^^ [I,!!, @| multiplied 
by a factor K{t) between and 1, namely the transmission coefficient, as we have got here 



K{t) 



detr 



m22 



cm 

^22 



detm, 



detm, 



(12) 



which describes the probability of a particle successfully 
escaped from the metastable well to recross the barrier. 
This expression for the fractional reactive index leads im- 
mediately to Kramers formula for the rate constant. 

In order to individually inspect the dynamical correc- 
tions of K{t) to the TST rate. Fig. Ogives the transient 
expression of it for various strength of system friction. 
Apparently seen from it, the value of K{t) keeps 1.0 when 
t — 0, showing no barrier recrossing behavior at the be- 
ginning. But as i — > 00, K(t) stabilizes to a invariable 
constant between and 1, namely the stationary trans- 
mission coefficient Kst, representing the probability of a 
particle already escaped from the metastable potential 
well to recross the barrier. Mathematically, it creates a 
cutoff factor to the TST rate. 

Meanwhile, by investigating the stationary transmis- 
sion coefficient at various strength of system friction we 
can find: in most cases of non-Ohmic damping, K{t) sta- 
bilizes to a smaller Kgt than the Ohmic case {S = 1.0) 
except for some special one such as 6 = 1.2. This implies 
the barrier recrossing of non-Ohmic system is generally 
stronger than the Ohmic case. The long memory aspect 



of friction is most possible to be at the origin of this 
behavior. 

For further understanding of this subject, the station- 
ary transmission coefficient Hgt is plotted in Fig. [5] as 
a function of S for various system frictions. Again it is 
found to change non-monotonically with 5. This is be- 
cause, the friction of non-Ohmic system has a long mem- 
ory in time. If the friction strength is weak (i.e. 6 —> 0), 
the motion of a particle depending greatly on the ini- 
tial conditions. In the opposite aspect, strong friction 
strength (i.e. S 2.0) causes great viscous force to the 
particle. In both cases, it is unfavorable for the parti- 
cle to pass the saddle point of the PES. This results in a 
smaller passing probability and a transmission coefficient 
than the normally Ohmic case. 

Furthermore, comparing with the results of ID non- 
Ohmic systems 1 13. 24], the 2D non-Ohmic system is 
found to have a bigger stationary transmission coefficient 
in the long time limit which means a weaker barrier re- 
crossing effect. Particle diffusing in the 2D PES under 
non-Ohmic friction seldom returns after passing the sad- 
dle point of the potential. Macroscopically it results in a 
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FIG. 4: Transient expression of K,{t) for various strength of 
system friction. Identical parameters are used as in Fig. [1] 
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FIG. 5: Stationary transmission coefficient function of 

S for various system frictions. Parameters used are identical 
to those in Fig (2] 

big net flux rate. 

IV. SUMMARY 

In conclusion, a two-dimensional barrier passage has 
been investigated, where the particle is subjected to a 
non-Ohmic memory friction with power exponent 5. The 
time-dependent passing probability and transmission co- 



efhcient are analytically obtained by using of the reaction 
flux method in a Langevin dynamic picture. In the long 
time limit, the passing probability, the optimal incident 
angle of the particle and the anomalous transmission co- 
efficient are all found to be non-monotonic functions of 
the power exponent of the friction S. The long memory 
aspect of friction is most possible to be at the origin of 
this behavior. The relatively bigger stationary transmis- 
sion coefficient reveals that there is less barrier recrossing 
in the two-dimensional system than the one-dimensional 
case. It is very difficult for the particle diffusing in the 
2D PES under non-Ohmic friction to return after passing 
the saddle point of the potential. Thus results in a big 
net flux rate. 
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APPENDIX. THE EXPRESSIONS OF {xi{t)) 

The quantities appear in the expression of {xi (t)) are 

P(s) = (mn/322[s]+m22/3ii[s]-2mi2/3i2[s])s3 

-|-(det/3ij[s] -I- miicj22 + '^22^^11 - 2mi2Wi2)s^ 

+ 0ii[s]uJ22 + (322[s]uJll - 2/3l2[s]t^l2)s 

+ {detmij)s'^ + detwy, 
Fi(s) = (mii/322 [s] + m-22/3ii[s] - 2mi2/3i2 [s])s^ 

-|-(det/3ij[s] -I- miiCi;22 - '^12^12)5 + (detmij)s^ 

+Pll[s]uj22 - /3l2[s]wi2, 
F2(s) = (mi2CJ22 - m-22t^l2)s + /3i2[s]lJ22 - /322[s]cJi2, 

Fsis) = (detmij)s^ + (toii/322[s] - mi2/3i2[s])s 
+miiUJ22 - "^12W12, 

¥4(3) = (mi2/322[s] - m22/3l2[s])s -I- 77112^22 - "^22C^12, 
F5(s) = m22S^ + /322[s]s + LJ22, 

Feis) = -mi2S^ - /3i2[s]s - t^i2- 

where /3ij[s] — ^ijS^^^ is the Laplace transform of the 
friction kernel function and jij — 7ija;J;^''sin^^((57r/2) is 
namely the effective friction. 
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